Infrared Spectroscopy of Quantum Crossbars 
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Infrared (IR) spectroscopy can be used as an important and effective tool for probing periodic 
networks of quantum wires or nanotubes (quantum crossbars, QCB) at finite frequencies far from the 
Luttinger liquid fixed point. Plasmon excitations in QCB may be involved in resonance diffraction 
of incident electromagnetic waves and in optical absorption in the IR part of the spectrum. Direct 
absorption of external electric field in QCB strongly depends on the direction of the wave vector q. 
This results in two types of ID — > 2D dimensional crossover with varying angle of an incident wave 
or its frequency. In the case of QCB interacting with semiconductor substrate, capacitive contact 
between them does not destroy the Luttinger liquid character of the long wave QCB excitations. 
However, the dielectric losses on a substrate surface are significantly changed due to appearance of 
additional Landau damping. The latter is initiated by diffraction processes on QCB superlattice 
and manifests itself as strong but narrow absorption peaks lying below the damping region of an 
isolated substrate. 
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I. INTRODUCTION 

Recent achievements in material science and technol- 
ogy have led to creation of an unprecedented variety of 
artificial structures that possess properties never encoun- 
tered in "natural" quantum objects. One of the most 
exciting developments in this field is fabrication of two- 
dimensional (2D) networks by means of self-assembling, 
etching, lithography and imprinting techniques 1,2 . An- 
other development is the construction of 2D molec- 
ular electronic circuits 3 where the network is formed 
by chemically assembled molecular chains. Such net- 
works have the geometry of crossbars, and bistable con- 
formations of molecular chain may be used as logical 
elements 4 . Especially remarkable is a recent experimen- 
tal proposal to fabricate 2D periodic grids from single- 
wall carbon nanotubes (SWCNT) suspended above a di- 
electric substrate 5 . The possibility of excitation of a 
SWCNT by external electric field together with its me- 
chanical flexibility makes such a grid formed by nan- 
otubes a good candidate for an element of random access 
memory for molecular computing. 

From a theoretical point of view, such double 2D grid 
i.e. two superimposed crossing arrays of parallel con- 
ducting quantum wires 6 " 8 or nanotubes 9 , represent a 
unique nano-object - quantum crossbars (QCB). Its spec- 
tral properties can not be treated in terms of purely ID or 
2D electron liquid theory. A constituent element of QCB 
(quantum wire or nanotube) possesses the Luttinger liq- 
uid (LL) like spectrum 10,11 . A single array of parallel 
quantum wires is still a LL-like system qualified as a 
sliding phase 9 provided only the electrostatic interaction 
between adjacent wires is taken into account. If an inter- 
wire tunneling is possible, the electronic spectrum of an 
array is that of 2D Fermi liquid (FL) 12,13 . 



Similar low-energy, long- wave properties are character- 
istic of QCB as well. Its phase diagram inherits some 
properties of a sliding phases in case when the wires 
and arrays are coupled only by capacitive interaction 9,14 . 
When inter-array electron tunneling is possible, say, 
in crosses, dimensional crossover from LL to 2D FL 
occurs 15,9,16 . If tunneling is suppressed and the two ar- 
rays are coupled only by electrostatic interaction in the 
crosses, the system possesses the LL zero energy fixed 
point, and a rich Bose-type excitation spectrum (plas- 
mon modes) arises at finite energies in 2D Brillouin zone 
(BZ) 14,19 . These QCB plasmons can be treated as a set of 
dipoles distributed within QCB constituents. In a single 
wire the density of the dipole momentum is proportional 
to the LL boson field 6(x) (x is the coordinate along the 
wire). 

Two sets of coupled ID dipoles form unique system 
which posesses the properties of ID and 2D liquid de- 
pending on the type of experimental probe. Some possi- 
bilities of observation of ID — > 2D crossover in transport 
measurements were discussed in Ref. 9. Here we con- 
sider various possibilities of direct observation of plas- 
mon spectra at finite frequencies and wave vectors by 
the methods of IR spectroscopy. 

In transport measurements, the geometrical factors 
regulate the crossover from anisotropic to isotropic resis- 
tivity of QCB: one may study the dc response for a field 
applied either parallel to one of the constituent arrays 
or in arbitrary direction. One may also study spatially 
nonuniform response by means of two probes inserted at 
different points of QCB and regulate the length scale i.e. 
the distance between two probes in comparison with the 
periods of the crossbar superlattice. These methods give 
information about nearest vicinity of LL fixed point at 
(g,w,T)->0. 

Unlike transport measurements, the methods of in- 
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frared spectroscopy provide an effective tool for inves- 
tigating the excitation spectrum in wide enough (q, uj) 
area well beyond the sliding phase region. We will show 
that the 1R spectroscopy allows scanning of 2D Brillouin 
zone in various directions and thereby elucidates dimen- 
sional crossover in the high symmetry points of the BZ. 

Several crossover effects such as appearance of non-zero 
transverse space correlators and periodic energy transfer 
between arrays ("Rabi oscillations") were discussed in 
our earlier publications 17 ' 18 . The direct manifestation 
of dimensional crossover is the response to an external 
ac electromagnetic field. To estimate this response one 
should note that the two main parameters characteris- 
ing the plasmon spectrum in QCB are the Fermi ve- 
locity v of electrons in a wire and the QCB period a 
(we assume both periods to be equal). These parame- 
ters define both the typical QCB plasmon wave numbers 
Q = IqI ~ Q = 2n/a and the typical plasmon frequen- 
cies uj <~ ojq = vQ. Choosing according to Refs.[ 11, 
5] v w 0.8 • 10 6 m/sec and a w 20 nm, one finds that 
characteristic plasmon frequencies lie in the far infrared 
region u ~ 10 14 sec -1 , while characteristic wave vectors 
are estimated as q ~ 10 6 cto _1 . 

In this paper we study high frequency properties of 
the simplest double square QCB (generalization to more 
complicated geometries is straightforward). We start 
from QCB interacting with an external infrared radia- 
tion. The plasmon velocity v is much smaller than the 
light velocity c and the light wave vector k is three or- 
ders of magnitude smaller than the characteristic plas- 
mon wave vector Q corresponding to the same frequency. 
Therefore, infrared radiaton incident directly on an iso- 
lated array, can not excite plasmons at all (it could excite 
plasmon with uj ^ 0). However in QCB geometry, each 
array serves as a diffraction lattice for its partner, giving 
rise to Umklapp processes of wave vectors nQ, n integer. 
As a result, excitation of plasmons in the BZ center q = 
with frequencies nvQ occurs. 

To excite QCB plasmons with q one may use an 
additional diffraction lattice (DL) with period A > a 
coplanar to the QCB. Here the diffraction field contains 
space harmonics with wave vectors 2itM/A, M integer, 
that enables one to scan plasmon spectrum within the 
BZ. Dimensional crossover manifests itself in the appear- 
ance of additional absorption lines when the wave vector 
of diffraction field is oriented along specific directions. In 
the general case one observes the single absorption lines 
forming two sets of equidistant series. Instead of that, in 
the main resonance direction (QCB diagonal) an equidis- 
tant series of split doublets can be observed. In the case 
of higher resonance direction, absorption lines form an 
alternating series of singlets and split doublets demon- 
strating new type of dimensional crossover related to the 
frequency change with direction fixed. 

We study also QCB interacting with a semiconducting 
substrate. This interaction may be rather strong because 



surface plasmon waves exist in the same frequency and 
wave vector area as plasmons in QCB. This problem is 
interesting from at least two points of view. First, tech- 
nological reasons may lead to creation of QCB embedded 
into thin layer of semiconductor medium. Second, the 
question of stability of LL regime in QCB against inter- 
action with 2D plasmons in substrate, is itself of great 
interest. It will be shown that in spite of the long range 
character of 2D screening, the substrate plasmons only 
renormalize the velocity of the QCB plasmons and do 
not destroy LL character of the QCB spectrum. At the 
same time, QCB-substrate resonance can be observed ex- 
perimentally in special directions. Moreover, at a special 
frequency, a triple resonance involving both QCB modes 
and substrate plasmons is possible. 

The structure of the paper is as follows. In section II, 
we briefly describe double square QCB and introduce the 
necessary definitions. Interaction of QCB with external 
field is studied in section III. In its first part we con- 
sider the case when the incident infrared radiation falls 
directly on the QCB (subsection III A). In the second 
part IIIB we study possible scanning of QCB spectrum 
with the help of an external DL. Dielectric properties of 
a combined system QCB-substrate are studied in the last 
section IV. In the Conclusion we summarize the results 
obtained. 



II. DOUBLE SQUARE QCB 

A square QCB is a 2D grid, formed by two periodically 
crossed perpendicular arrays of ID quantum wires or car- 
bon nanotubes. Arrays are labeled by indices j = 1, 2 and 
wires within the first (second) array are labeled by an in- 
teger index (ni). In experimentally realizable setups, 
QCB is a cross-structure of suspended single- wall carbon 
nanotubes lying in two parallel planes separated by an 
inter-plane distance d, placed on a dielectric or semicon- 
ducting substrate (see Fig.l). Nevertheless, some generic 
properties of QCB may be described under the assump- 
tion that QCB is a genuine 2D system. We choose coor- 
dinate system so that 1) the axes Xj and corresponding 
basic unit vectors are oriented along the j-th array; 2) 
the X3 axis is perpendicular to QCB plane; 3) x 3 coor- 
dinate is zero for the second array, — d for the first one, 
and —(d + D) for the substrate. The basic vectors of 
the reciprocal superlattice for a square QCB are Qei ^, 
Q = 2ir/a so that an arbitrary reciprocal superlattice 
vector m is a sum m = mi + m.2, where m.j = rajQej, 
(raj integer). The first BZ is a square \q\^\ < Q/2. 
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FIG. 1. QCB on a substrate, ej - basic vectors of the co- 
ordinate system. Vector ei (e 2 ) is oriented along the first 
(second) array. The inter-array distance is d and the distance 
between substrate and the first (lower) array is D. 

A single wire is characterized by its radius r , length L, 
and LL interaction parameter g. The minimal nanotube 
radius is r w 0.35 nm 20 , maximal nanotube length is 
L « 1 mm, and the LL parameter is estimated as g rts 
0.3 11 . In typical experimental setup 5 the characteristic 
lengths mentioned above have the following values 

d Rj 2 nm, L w 0.1 mm, 

so that the inequalities 

ro«d«a«L 

are satisfied. 

The Hamiltonian of QCB interacting with an external 
field is 



H = Hqcb + He- 



(1) 



The QCB interaction He with an external electric field 
E = (Ei, E 2 , E 3 )) is nothing but an energy of a set of 
QCB dipoles in this field 



Hp 



, n 2 a, —d)9i(x 1 ,n 2 a, -d)+ 

{ n 2 J 

^2 J dx 2 E 2 (n 1 a,x 2 ,0)9 2 (n 1 a,x 2 ,0), > 



where 9j is one of the two conventional canonically con- 
jugate boson fields tvj, 9j. 



The QCB Hamiltonian 

Hqcb = Hi + H 2 + H\ 2 . 



(2) 



consists of three terms. The first of them describes LL 
in the first array 

L/2 

H\ = -y y~] / dx 1 lgTrl(x 1 ,n 2 a, -d) 

1 2 ) 

+ - (d Xl 9i (xi,n 2 a, -d)) >. 

The Hamiltonian H 2 is obtained from Hi after permuta- 
tion 1 <-> 2 and replacement — d — > in the arguments of 
the fields. 



The inter-array interaction is described by the last 
term in Eq.(2) 



#12 = Vo ^2 / dxidx 2 ( f 



2e 2 



It results from a short-range contact capacitive coupling 
in the crosses of the bars. The dimensionless envelope 
function (introduced phenomenologically) C( r j) describes 
re-distribution of a charge in a tube j induced by the 
interaction with tube i. This function is of order unity 
for |t| ~ 1 and vanishes outside this region so that the 
dimensionless integral 



/ 



C(r)e 



ikroT 



dr ~ 1 



xd xi e 1 (x 1 ,n 2 a,-d)d X2 6 2 (n 1 a,x 2 ,0), V = —. (3) 

d 



is of order unity for all |fe| smaller than a certain ultra- 
violet cutoff. For simplicity, we put it equal to unity in 
what follows. 

Before turning to investigation of QCB interaction 
with an external field, we briefly describe the spectral 
properties of the QCB itself 18 . The QCB Hamiltonian 
(2) is a quadratic form in terms of the field operators, 
so it can be diagonalized exactly. Such a procedure is 
rather cumbersome. However, due to the separability of 
the interaction (3) the spectrum can be described ana- 
lytically. Moreover, being small, this interaction grossly 
conserves the unperturbed ID systematics of levels and 
states, at least in the low energy region corresponding to 
the first few energy bands. This means that perturbed 
eigenstates could be described in terms of the same quan- 
tum numbers (array number, band number and quasimo- 
mentum) as the unperturbed eigenstates of an "empty" 
lattice. Such a description fails in two specific regions of 
a reciprocal space k = fciei + fe^- The first of them 
is the vicinity of high symmetry lines kj = nQ/2 with n 
integer (the lines with n = ±1 include BZ boundaries). 
Around these lines, the interband mixing is significant 
and two modes from the adjacent bands are degenerate. 
The second region is the vicinity of the resonant lines 
k\ ± k 2 = nQ where the eigenfrequencies of unperturbed 
plasmons from the same band s propagating along two 
arrays coincide u>i s (qi) = Lu 2s (q2)- Around resonant lines 
inter- array mixing is significant and at these lines two 
modes corresponding to different arrays are degenerate. 

Inter-array interaction introduces real two dimension- 
ality into the problem, lifts the degeneracy of bare modes 
and splits degenerate frequencies. Therefore the simplest 
way to probe 2D nature of QCB is to observe this split- 
ting in the corresponding spectral region. Possible ways 
of such an observation are discussed in the next two sec- 
tions. 



III. INFRARED LIGHT ABSORPTION BY QCB 

A. Long Wave Absorption 

In the case of a dielectric substrate transparent in the 
infrared region, one can treat QCB as an isolated grid 
(without substrate) interacting directly with the incident 
radiation. Consider the simplest geometry (see Fig. 2 for 



3 



details) where an external wave falls normally onto QCB 
plane, and its electrical field 

E = E ei cos (kr — tot) 

is parallel to the lower (first) array. In this geometry 
the field E is longitudinal for array 1 and transverse for 
the array 2. The cigenfrequencies of transverse modes 
in array 2 substantially exceed the IR frequency of the 
incident wave and even the standard LL ultraviolet cut- 
off frequency. Thus, the incident wave can be treated 
as a static polarization field for this array, and the fac- 
tor costot can be omitted. Then, the polarization waves 
in array 2 form a longitudinal diffraction field for array 
1 with quasi wave vectors nQ (n integer). Further, the 
characteristic order of magnitude Q of a QCB plasmon 
wave vector is much larger than the wave vector k of 
the incident light, and we put the latter equal to zero 
from the very beginning. Then, the light wavelength is 
much larger than a nanotube diameter and the geometri- 
cal shadow effect can be neglected. As a result the total 
field which affects array 1 consists of an external field and 
a diffraction field produced by a static charge induced in 
array 2. 




FIG. 2. The incident field orientation with respect to QCB. 
The axes x\ and X2 are directed along the corresponding ar- 
rays, and d is the inter-array vertical distance (along the xs 
axis). 

To calculate diffraction field consider first the field E° 
produced by the quantum wire of array 2 which is lo- 
cated at x\ — x 3 = and labeled by n\ = 0. The large 
distance between the wire under consideration and its 
neighbour partners from the same array allows us to ne- 
glect the influence of the charges induced on them. The 
static potential on the surface of the wire includes ex- 
ternal potential of an incident field and the potential $° 
of the charge induced on the wire. On the other hand 
this static potential should be equal to a constant which 
we choose to be zero. In cylindrical coordinates r, $,X2, 
x\ = rcosi9,x 3 = rsin$, this condition reads 

$°(R ,'&,x 2 ) =E r cos#. (4) 

Outside the wire, induced potential $° satisfies Laplace 
equation A<E>° = 0. Solving this equation with boundary 



condition (4) wc obtain the static part of the induced 
potential 

$°M,x 2 ) = ^^costf 
r 

and the corresponding static part of the induced field 
along x\ direction 



E°{ Xl ,x 3 ) = -Eq 



r o (4 - x i) 

(xl + xlf ■ 



The first component of the difraction field is the sum of 
the fields induced by all wires of the upper array, 

Ei(xi; t) — cos cot Ei (x\ — n^a, —d) 

711 

tp ^ 4 (d 2 - (xi - ma) 2 ) 

= -E a cosLoty — tt~ . 5 

{d 2 + {x 1 -n 1 aYf 

This field is a periodic function of x\ with period a. 
Therefore, its Fourier expansion contains only wave vec- 
tors ki n = nQ (n is the order of diffraction) . This means 
that only frequencies w„ = nvQ can be excited. In this 
case it is more convenient to expand the field over Bloch 
eigenfunctions of an "empty" wire 18 . These functions are 
labeled by quasimomentum qi, \qi\ < Q/2, and the band 
number s. The expansion includes only qi = compo- 
nents and has the form 

Ex{xi;t) = cosLot^2E [s/2 ]U s (xi), 



where 



,(x) =cxp (iQx[s/2](-l) 



is the qi = + Bloch amplitude u sqi (x) within the s-th 
band, [. . .] is the entire part symbol, and 



E n = -E ^nQde- nQd . 
ad 



(6) 



The excited eigenfrequency ut n — ui[ s /2] belongs simulta- 
neously to the top of the lower even band with number 
s = 2n and to the bottom of the upper odd band with 
number s = 2n + 1 (this is the result of E(x) parity). 
Incident field corresponds to n — and we do not take it 
into account. 

Turning to q, s representation with the help of the ex- 
pansion 

6i(xi,n 2 a) = ^J2^e i(9ixi+q2n2a) 



x cxp (iQx 1 [s/2] (-l) s ^ign^) 



(7) 



and similarly for 82 and ni.2, one easily sees that only the 
q = components are involved in interaction with the 
incident radiation. The corresponding part of the Hamil- 
tonian (1) has the form ( zero quasimomentum index is 
omitted): 
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where 



hva 



0.007 



is the dimcnsionless inter-array interaction. Consider the 
initial frequency oj close to ut n . In a resonant approxima- 
tion, only four equations of motion for the "coordinate" 
operators 6 S with s — 2n, 2n + 1 are relevant 

0l,2n + ^l8l,2n + fan (#2,2n - 2 ,2n+\) = Lf n COS U)t, 

6l,2n+l + <^n#l,2ra+l — fan (&2,2n — #2,2n+l) = Lf n COSUjt, 

#2,2n + Un&2,2n + fan (#l,2n ~ #l,2n+l) = 0, 

#2,2n+l + <^#2,2ra+l - fan (#l,2n - 01,2ra+l) = 0, (8) 

where 

_ V2vge 

The homogeneous part of this system defines four eigen- 
frcquencies 

^99 = W «3 = 
U uu ~ W„(l - 0), 

U gu ~ W n (l + </>)• 

The corresponding eigenvectors are symmetrized combi- 
nations of the four operators which enter Eq.(8). They 
have a fixed parity with respect to permutation of arrays 
(the first index) and neighboring bands (the second in- 
dex) . Only two modes (even with respect to band index) 



®99 ~ o (^b 2 ™ + 01,2n+l + #2,2n + #2,2n+l) > 



^"9 — o (^!> 2 ™ + #1,271+1 — #2,2n — ^2,2n+l) 



interact with an external field. Therefore only the unper- 
turbed frequency u> n = u> gg = uj ug will be absorbed. The 
two equations of motion for the operators 9 ggyU g have the 
same form 

9 a + 2j6 a + u n 6 a = Lf n cos ujt, 

where a = gg, ug. Employing standard procedure in the 
vicinity of the resonance \u) — u) n \ <C 0J n immediately 
yields the relative absorption of Lorcntz type 



where 

is the energy of light that falls on the QCB per unit time. 

Due to the exponential term in the r.h.s of Eq.(6), E n 
decreases fast with n and only the first few terms con- 
tribute to absorption. The characteristic dimensionless 
scale of the induced field r^Kad) for typical values of 
QCB parameters equals 0.004. We tabulate below the 
lowest dimensionless Fourier components of the induced 
field. 



n 


1 


2 


3 


4 


5 


ad ti n 


1.05306 


1.12359 


0.89914 


0.63957 


0.42650 



The results show that one can hope to probe at least 
the first five spectral lines corresponding to u n with 
n = 1,2, ... ,5. 

The width of the absorption line (9) is governed by an 
attenuation coefficient 7. It was introduced phenomeno- 
logically but one can (at least qualitatively) estimate its 
value. The attenuation is caused by decay of plasmon 
into phonons. The one phonon decay of the plasmon 
with wave number k and frequency lu = v\k\ into a single 
phonon with the same u) and k occurs in a single point 
in ID and does not yield finite attenuation at all. Mul- 
tiphonon decay is weak because of the small anharmonic 
coupling within the wire. As a result, the form of the 
absorption lines should be determined mainly by the in- 
strumental linewidth. 



B. Scanning of the QCB Spectrum within BZ 

Within a geometry considered in the previous subsec- 
tion, one can probe plasmon spectrum only at the BZ 
center. To study plasmons with nonzero wave vectors 
one should add to the system an external diffraction lat- 
tice namely a periodic array of metallic stripes parallel 
to the Y axis (see Fig. 3). The DL plane Z = is parallel 
to the QCB planes Z = —D for the upper second array 
and Z = —(D + d) for the lower first array (the Z axis 
is parallel to the X3 axis). The distance D between DL 
and second array is of the same order that the inter-array 
distance d = 2 nm. The angle between DL wires and the 
second array is ip (0 < ip < tt/2). To get a wave number 
if of a diffraction field much smaller than Q one needs 
a DL with a period A much larger than the QCB pe- 
riod a. In the following numerical estimations we choose 
A w 200 nm. 
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[nQde~ nQd }\ (9) 
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FIG. 3. QCB and DL. The X, Y axes are oriented along 
the DL stripes and the wave vector K of the diffraction field 
respectively. The DL (QCB) period is A (a). 

Consider an incident field with electric vector E = 
Eq&x cos(kr — cut) oriented along the X axis (perpen- 
dicular to the DL wires). The radius Ro of a DL wire is 
assumed to be not much larger than the nanotube radius 
tq. In this case light scattering on the DL is similar to 
that considered in subsection III A) . Then the diffraction 
field is concentrated along the X direction and has the 
form (compare with Eq.(5)) 



E X (X, Z,t) -- 

Rl{z 2 



— Eq cos u>t x 
(X — NA) 2 ) 



N 



(Z 2 + (X - NA) 2 ' 
The Fourier transform of the diffraction field is 
E X (K,Z) 

where 

K(M) = Ke x = (K!,K 2 ) 



2ttM 



-\kz\ 



(sin ip, cos if) 



with positive integer M. This means that all the points K 
lie on the same ray oriented along the positive direction 
of the X axis. The vector K(M) for a fixed M can be 
uniquely represented as a sum of quasimomenta lying in 
the first BZ and two reciprocal lattice vectors 

K(M) = q(M) + mi(M) + m 2 (M). 

The field components 



E 1VL = E x {K,D + d)sm<p, 
E 2 k = E x (K, D) cos ip 



(10) 



parallel to the quantum wires, can excite plasmons and 
contribute to the absorption process. 

The Hamiltonian (1) of QCB interacting with an ex- 
ternal field in Fourier 

■ Z hr ( ^-- "-•■»}• i = 1 ' 2 ' 



Hl2 = ^Yl $ q+m#i, q+m Aq+m 2 
m,q 

E[/i,k (oi, K + e\ K ) 

K 

+/2,K ^2,K + ^2,K 



(12) 



He — 



TiL 



+ 



where 



$k = 0£fci6fc 2 ! $fc=w fe signfe, uj k =v\k\, 
V2vge 

• /: ' K i^r'-' K - m 1111 1 m - 



(13) 



In this subsection we are interested not in the form of 
the absorption line but only in the resonant frequencies. 
Therefore, we did not introduce any phenomenological 
attenuation. The equations of motion for boson fields 
have the form 

£*l,q+mi + w ? 1 +m 1 Q^l,q+mi + 

,q+m 2 — 

m 2 M,m 2 
,2 



^2,q+m 2 + ^g 2+m2 Q#2, q +m 2 + 
+ ^ ^q+m^Lq+mi =L ^ /2,K(5K,q+m. 



(14) 



M,mi 



Only the first few terms in the sum over K in the r.h.s. 
of Eq. (14) really excite the QCB plasmons. Indeed, 
the diffraction field (7) is proportional to the same di- 
mensionless function of the type te^ t (t = \KZj\) as in 
the previous subsection (see Eq.(6)). This function has 
its maximum at t = 1 and differs significantly from zero 
for 0.2 < t < 2.7. For a = 20 nm, D = 2 nm, it is 
of order unity within the interval 0.18Q < \K\ < 2.13Q 
for the first array (Z\ = D + d), and within the interval 
0.36Q < \K\ < 4.26Q for the second array (Z 1 = D). 
This means that one can excite the modes of four lower 
bands (K < 2Q) of the first array and the modes of 
eighth lower bands (K < AQ) of the second array. 

According to Eqs. (10) the field iyK(M) is coupled 
with plasmons of wave vectors q + u\j — q(M) + mj (M) 
within the j-th array. The nature of the excited plas- 
mons as well as their frequencies depend on the direction 
of the vector K(M). For simplicity we restrict ourselves 
by acute angles < if < ir/2 describing orientation of 
both DL and vector K(M). There are four kinds of di- 
mensional crossover depending on specific directions in 
the BZ. Each type of crossover is characterized by its own 
set of absorption lines. The first one takes place in a com- 
mon case when K(M) for any M never reaches neither a 
resonant direction nor the BZ boundary. The second case 
corresponds to the bisectorial direction ip = 7r/4 where 
the main resonant condition lu(Ki) = lu(K2) is fulfilled. 
The third group of directions is determined by another 
resonant condition oj(Ki) = u(nQ =F K 2 ). Finally, the 
fourth group is formed by directions intersecting with 
the BZ boundaries for some values of M. In what follows 
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we consider these four cases separately. 

1. In the general case, the points K(M) for all M 
are far from the BZ diagonals and boundaries. There- 
fore each of them corresponds to two plasmons mostly 
propagating along the j-th array, j — 1,2, with unper- 
turbed frequencies oJk^m) = v Kj{M). The inter-array 
interaction slightly renormalizes the eigenfrcqucncies 
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Thus, increasing the frequency of an incident light one 
observes a set of single absorption lines that consists of 
two almost equidistant subsets with frequencies corre- 
sponding to excitation of plasmons in the first or second 
arrays. The distances between adjacent lines within each 
subset are 

Au>i — vAKi — 2irv sin <p/A, 
Au 2 = vAK 2 = 2nvcosLp/A, 

and their ratio depends on the DL orientation ip only 



Alu 2 



= tan ip. 



2. In the resonant case <p = 7r/4, the relation K\{M) = 
K 2 (M) is satisfied for all M . Therefore modes propagat- 
ing along the two arrays are always degenerate. Inter- 
array interaction lifts the degeneracy. Indeed, in the res- 
onant approximation, the coupled equations of motion 
for the field operators read 

#1K + W^-^ik + 4' u, k 1 2 k. — /lK, 
#2K + w|- 1 02K + fp^K^lK = ./2K- 

After symmetrization 9 g ^ u — (#ik ± ^kVv^, they have 
the same form 



9a + ^ a a — fa, 



where 



are the renormalized frequencies, a — g,u, and f 9iU = 
(/ik ± f 2 n)/\/2- The amplitudes f g . u are of the same 
order of magnitude because the distances D and d are 
different but have the same order of magnitude. As a 
result, increasing the frequency of an incident light one 
observes an equidistant set of absorption doublets with 
distance ny/2v/A between adjacent doublets. 

3. Consider now the directions tp determined by the 
equation 



nA 



V2M a 



where n and Mq are mutually prime integers. For this di- 
rection, two components of the first M — 1 points K(M) 
do not satisfy any resonant condition while the M -th 
one does 



K 1 (M )±K 2 (M ) = nQ. 



(15) 



With increasing M this situation is reproduced periodi- 
cally so that all points K(pMo) with p integer satisfy a 
similar condition with pn standing instead of n, while all 
intermediate points are out of resonance. 

In the zero approximation with respect to the inter- 
array interaction we expect to observe two set of ab- 
sorption lines with frequencies pujj = vKj(pMo), j — 
1,2, corresponding to excitation of plasmons within the 
pnij(Mo)-th band of the j-th array. The ratio of the 
frequencies ojj is defined by DL orientation 



uo 2 



= tan<p. 



However, due to the resonance condition (15), plasmon 
in the first array with wave vector Ki(pMo) and fre- 
quency uji = vKi(pM ) is coupled with plasmon in the 
second array with the same frequency and wave vec- 
tor K' 2 = T{ n pQ — Ki(pMoj) (inter- array degeneracy). 
Similarly, plasmon in the second array with wave vector 
K 2 (p Mo) and frequency uj 2 = vK 2 (pMo) is coupled with 
plasmon in the first array with the same frequency and 
wave vector K[ = npQ =F K 2 (pMo). This degeneracy of 
two modes corresponding to the same band but to dif- 
ferent arrays is lifted by the inter-array interaction. As 
a result one has two sets of doublets instead of two sets 
of single lines. 

Thus, for such orientation of DL, increasing frequency 
of an incident wave one should observe two equidistant 
sets of single absorption lines with two sets of equidistant 
doublets built in these series 

1 

T 
i 

2* 

In the case n = 1 the lower doublet lies in the first en- 
ergy band, whereas the upper one lies in the second band. 
For A/ a = 10 (that corresponds to the realistic values of 
the parameters a = 20 nm and A — 200 nm) the lowest 
doublet (p = 1) will be observed for example for integer 
M = 8, at the angle ip(8) « 17°, around frequencies 
wi(8) = 0.76i;Q, uj 2 (8) = 0.24vQ. 

4. It seems that a similar behavior will be manifested 
in the case when the points K p m lie at one of the BZ 
boundaries, i.e. satisfy the relation 



Wik — ^>K 1 (pM ) ^1 ± 
W 2 K = ^K 2 (pM ) 1 ± 



npQ 
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with some specific values j, Mj and n. Such situation 
is realized at specific angles that depend on the integers 
j,n,Mj. In the vicinity of the points K(pMj) two fre- 
quencies corresponding to the unperturbed modes of the 
j-th array from the np-th and (np + l)-th bands coin- 
cide. This is the case of inter-band degeneracy that is 
also lifted by inter-array interaction. Due to the square 
symmetry (invariance with respect to Xj — > — Xj inver- 
sion), only one of the two components with frequency 
cj = v\Kj(pMj)\ may be excited by a diffraction field. 
Therefore, this case does not differ from the case 1 con- 
sidered above and two sets of equidistant single lines can 
be observed. 

We emphasize that studying absorption of light by 
QCB one can expose, beyond the standard 18 dimensional 
crossover with respect to an angle (direction), also oc- 
currence of a new type of crossover with an external fre- 
quency as a control parameter. This occurs for special 
directions of type 3 where, with increasing frequency, the 
set of single lines is periodically intermitted by doublets. 



IV. QCB ON A SEMICONDUCTING SUBSTRATE 

A. Hamiltonian 

An alternative method of examination of the plasmon 
spectrum in QCB is the study of its interaction with 
the substrate, whose electrodynamic properties are well 
known. An excellent material for this purpose is a III-V 
semiconductor. Surface plasmons in semiconductor slab 
can be excited in a controllable way (see Refs. 21,22 and 
literature therein), and these plasmons can be in reso- 
nance with Bose excitations in QCB. 

In this section we consider dielectric properties of QCB 
placed on a semiconducting substrate (Fig.l). Any sur- 
face wave excited in the substrate is coupled with QCB- 
plasmon modes due to the substrate-QCB interaction. 
Such an interaction gives an opportunity to probe QCB 
spectrum via IR absorption of the substrate (see below). 
Assuming the distance D between the first array and the 
substrate to be much smaller than the distance d be- 
tween arrays, one can keep only the interaction between 
the substrate and the first array. In this case the system 
is described by the Hamiltonian 

H = H a + H mt . (16) 

Its unperturbed part 

Ho = Hk + H i + H 2 

contains the Hamiltonians (11) of the two arrays and the 
kinetic energy of the substrate electrons with effective 

mass m and quadratic dispersion law = — — (we 

2m 

omit the irrelevant spin variables) 



Hk — £q+mCq +m C q+m . 



Here as in the previous section, a quasi wave vector q be- 
longs to the first BZ and m is a reciprocal lattice vector. 
The interaction part of the Hamiltonian (16) 



Hint — He + H s \ + H 
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includes inter-array interaction (12), Coulomb interac- 
tion within the substrate 



H c 



q+mPq+ m Pq+m i 



qm 



2^e 2 
k 



k = q + m 



and an interaction H s i between the substrate and the 
first array. The latter interaction is a capacitive cou- 
pling between charge fluctuations in the substrate and 
collective modes in the nearest array. In coordinate rep- 
resentation it is written as 

H s i = ^ / dxidx[dx 2 W(xi — x[,X2 — n 2 a) x 

xp(x 1 ,x 2 )d x ' i 6 1 (x' 1 ,n 2 a). 

Here p(r) is the density operator for the substrate elec- 
trons, \/2d x ' i 6i(x' 1 , n 2 a) is the density operator of the 
first array, and W(r — r') is Coulomb interaction between 
the substrate and array 



W(r) 



V2e 2 C P 
^f\r\ 2 + D2 



In momentum representation the interaction between the 
substrate and the array has the form: 



H s i — 



(17) 



mq 



where 



Wk = ik 1 J^- [ d Xl dx 2 W(r)e lkr 
V ha J 

is proportional to the Fourier component of IT(r). 



B. Dielectric Function 

The high frequency properties of the system at zero 
temperature are determined by the zeroes of the dielec- 
tric function 



1 



e(k,w) 



i + [/ k n(k,w). 



(18) 



Here 
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n(k,w) 



dte v 



(19) 



is the polarization of the substrate interacting with QCB, 
Pk(t) = e lHt / h pke~ lHt / R is the density of the substrate 
electrons in the Heisenberg representation, and averag- 
ing is performed over the ground state of the Hamilto- 
nian (16). 

Dielectric properties of the substrate per se within the 
RPA approach are described by the 2D version of the 
Lindhard formula 

e s (k,w) =l-[/ k n (k,o;), 

n , ,\ - 1 6<y - £F Z £k ') ~ e j£L Z £k + k ') (on\ 

no(k,«)-^^ ^-(^-^o + io • (20) 

The substrate polarization II s (k, u>) is defined by the same 
Eq.(19) where the averaging in the r.h.s. is performed 
with respect to the substrate Hamiltonian H s = Hk + 
He- It can be written as 

(II s (k, ^)) -1 = (n (k, ^)) -1 - U k , (21) 

that results in the expression for the dielectric function 

1 



e s (k,w) 
similar to Eq.(18). 



l + U k U s (k,uj) 



(22) 



The active branches of substrate excitations are the 
surface density fluctuations which consist of 2D electron- 
hole pairs and surface plasmon mode with dispersion law 



w s (k) = v F kJl + 



2kr B 



(23) 



h 2 

where r B = 7 (see e.g. Rcf. 23). The RPA spectrum 

of surface excitations is shown in Fig. 4. 

For GaAs, the substrate parameters are to = 0.068mo, 
too = 9.1 • 10 -28 g is the mass of a free electron, 
vf = 8.2 • 10 6 cm/sec, r B — 0.78 nm, ll>o = vf /r B = 
1.05 x 10 14 sec -1 . For k < k* rts O.lr^ 1 , the plasmon 
frequency lies above the continuum spectrum of electron- 
hole pairs and the substrate plasmons are stable. Besides, 
one may easily satisfy the resonance condition for the col- 
lective plasmon mode near the stability threshold k ~ k* 
and QCB excitations with frequency ~ 10 14 sec -1 . For 
large enough k > k* the plasmon dispersion curve lies 
within the quasi-continuum spectrum and plasmons be- 
come unstable with respect to decay into electron-hole 
pairs (Landau damping of the substrate plasmons). 



w/wo 
0.25 




0.05 0.1 0.15 0.2 
FIG. 4. Dispersion of the substrate plasmons (upper line) 

and quasi-continuum spectrum of electron-hole excitations, 

(dashed area). Frequency is measured in loo = vf/th units. 

To obtain the dielectric function of the whole system 
substrate plus QCB (18), one should take into account 
that Umklapp processes stimulated by interaction be- 
tween the substrate and the first array (17) as well as 
the interaction between arrays (12), produce modes with 
wave vectors q+m with various reciprocal lattice vectors 
m. This necessarily leads to appearance of non-diagonal 
components of the polarization operator 



n(q+m,q+m';cj) 



00 

/ 



dte 1 ' 



Pq+m(t)pq +m /(0) 



In what follows we always consider a fixed frequency to 
and a fixed wave vector q in the BZ. Henceforth the vari- 
ables q and to are omitted below for simplicity. In the 

framework of RPA approach, II(m, m') = "HMJo satis- 
fies the following Dyson-type equation 



in in 

<QTTTTTD> 



(24) 



Here the first term C Z> = H s (m) in the right hand side 

is the substrate polarization (21) of the isolated substrate 

itself, the empty point = W m = VF q + m is a simple ver- 
tex which describes substrate - (first) array interaction, 
and the thick line 



oc 

"j j m ' _ _l_ f 

h J 



dte 1 ' 



0; q +rn 3 (*),pt +m ,(O)l) (25) 



is the correlation function of the j-th array mode and the 
substrate plasmon. 

The Dyson equation (24) should be completed by two 
equations for the correlation functions (25) (j = 1,2) 




(26) 
(27) 



J m, 



Here the dashed line ■ — D®(rrij) (j = 1,2) is the 

bare correlation function of the j-th array modes 
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D°J mj ) 



OC 

-~f 

vg J 



dte 



7jq+m 



(*).^(0)]) Q 



uj 2 — v 2 (qj + raj) 2 



and two triangles A = y // h<j)/vg^ qi+mi Q and V = 
^h(j)/vg£ tq2 + rn2 Q form the second vertex describing the 
separable inter-array interaction (13). 

Solving the system of equations (24), (26) and (27) one 
obtaines the diagonal element n(m) = IT(m, m) of the 
polarization operator 



p(m)]- 1 - [iMm)]- 1 - |^ m | 2 D(m) 



(28) 



The second term in the right-hand side of this equation 
describes renormalization of this vertex by interaction 
between the substrate and QCB. The factor -D(m) is a 
renormalized correlation function of modes of the first 
array 

[^(m)]- 1 = [D?(mi)] _1 - Mm) + p(mi)). (29) 

Here w(m) describes the effective interaction between 
the first array and substrate 



w(m) =F(mi)- |^ m | 2 n s (m), 



F( mi ) 



m 2 



|W m | 2 n s (m) 



(30) 



and <p(mi) is the effective interaction between the two 
arrays 



<p(mi) 

u mj = v(qj +mjQ), 



m 2 



(31) 



renormalized by Coulomb interaction of array modes 
with substrate plasmons, 



4V 



J^ mi WK)) 



(32) 



Equations (28) - (32) together with the definition (18) 
solve the problem of dielectric properties of the combined 
system QCB-substrate within RPA. 

The whole {lu, k} space can be divided into two parts. 
Within the first subspace the dielectric function (18) is a 
real function. Its zeros define the spectrum of collective 
modes. Within the second subspace the dielectric func- 
tion has non-zero imaginary part. The latter is related 
to regions where these excitations are unstable. These 
two subspaces will be considered separately in the next 
two subsections. 



C. Excitations 

The spectrum of collective excitations in QCB- 
substrate system is determined by zeroes of the dielec- 
tric function e(q, u>) = 0. In the low frequency limit 
ui <C vQ/2 taking into account only contribution of the 
lowest QCB band (i.e. omitting the terms with m' ^ 0), 
we obtain the dispersion equation 



1(9)) 



(a; 2 -(l-a)A?) ( W 2 - (1 - a)« 2 «§) 



'12 



TT z n 

±2 a 



2 2\ 



(33) 



where u) s (q) is the substrate plasmon frequency (23). Ne- 
glecting both inter- array (oc 4> 2 ) and substrate- first array 
(oc W^) interactions, we are left with the substrate plas- 
mon mode and the modes of two arrays. Asymmetry of 
the dispersion relation (33) with respect to the two QCB 
arrays reflects the fact that the substrate interacts only 
with the first array. 

The key question is the robustness of the QCB spec- 
trum against interaction with 2D substrate excitations. 
To answer this question wc note that for wave vector 
q — > both interactions vanish 



,2 ,2 



I Wo 







(34) 



for all directions of q except the first array direction x\ 
(the latter direction is singular in all similar II? — 2D 
problems). Therefore in the long wave limit q -C Q the 
interaction just renormalizes the bare dispersion laws of 
the arrays 

= 1 - a - dljOtq, Ctq = 



2 n 2 



2nv 2 e 
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conserving its LL linearity. This result verifies stability 
of QCB plasmons with respect to substrate-QCB inter- 
action. 

The dispersion relation (33) shows that the QCB- 
substrate interaction results in appearance of two more 
resonant directions in the BZ. In addition to the diagonal 
line OJ (Fig. 5) that is the resonance line for inter-array 
interaction, two new lines arise, which correspond to reso- 
nant interaction of the substrate with the first or the sec- 
ond array. The resonance condition lo s (q) = vqj = uj (q) 
is fulfilled along the line LN for j = 1 and along the 
line KM for j — 2. The points J, M on these lines are 
the points of onset of Landau damping. The region of 
Landau damping exists outside the curve WJUMR and 
consists of two subregions. One of them exists in an iso- 
lated substrate outside the arc PNUMR. Another region 
is secluded within the close curve WJUNPW and ap- 
pears exclusively due to QCB-substrate interaction (see 
subsection IV D for details) . 
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0.1 0.2 0.3 0.4 0.5^ /Q 

FIG. 5. The first quarter of QCB BZ (a = 20 nm). QCB 

plasmons are in resonance with surface plasmons along the 
lines OM and ON. The surface plasmons are stable within 
the internal part of BZ bounded by the curve WJU and the 
arc UMR. 

Far from the symmetric point B, the two renormalized 
dispersion laws along LN are 



w 4 



= v 2 q\ (1 ± V^) • 



(35) 



They describe the splitting of two coupled substrate-first 
array modes. The third mode is slightly renormalized 
plasmon propagating along the second array. Along the 
resonance line KM one of the dispersion laws conserves 
its bare form 

lu = vq 2 

while another one transforms into 



„ 2 = 1 - (a q + a) - 2 



(36) 



These two modes are coupled substrate plus second array 
modes. The third mode represents slightly renormalized 
plasmon propagating along the first array. Along the 
diagonal OJ, the two array modes are in resonance vq\ = 
vq 2 = w(q). Here we get two coupled QCB plasmons with 
dispersion laws 



oj± = vq(l ± s/a) 



(37) 



and substrate plasmon with slightly renormalized fre- 
quency. At the symmetric point B, qis — q2B ~ 0.1Q, all 
three modes are in resonance w s (gs) = vqiB = vq 2 B = 
iv b- Interaction lifts the degeneracy and renormalizes the 
frequencies 



B ' 



| (l± V^5 



(38) 



In Figs. 6, 7 the dispersion curves calculated along the 
circle arcs with two different radii are shown. Solid lines 
show the dispersion curves of interacting system whereas 
the dashed lines describe those of noninteracting sub- 
strate and QCB. Here the azimuthal angle ip is intro- 
duced as in the previous section. 



0.2 0.4 0.6 0.8 1 
FIG. 6. Frequencies of substrate plasmons and modes of 

QCB with interaction between them (solid lines) and without 

the interaction (dashed lines) as functions of the azimuthal 

angle ip for q = 0.3Q. 

We start with a general case (Fig. 6) that represents 
the dispersion curves along the DEFGH line for q = 
0.3Q. Points D, E, F, G, H in Fig. 5 correspond to values 
2ip/ir = 0, 0.28, 0.5, 0.72, 1. Far from the resonance 
lines (i.e. not too close to the points E, F and G) the dis- 
persion curves describe slightly renormalized frequencies 
of the modes propagating within the substrate and along 
two arrays. At the point E in Fig. 5, the modes propa- 
gating in the substrate and along the first array are de- 
generate. The substrate-QCB interaction splits this de- 
generacy and the renormalized frequencies are described 
by Eq.(35). At the point F, the modes propagating along 
the first and second arrays are degenerate and the fre- 
quencies arc described by Eq.(37). At the point G, the 
modes propagating in the substrate and along the second 
array are degenerate and the frequencies are described 
by Eq.(36). For any fixed angle ip three absorption lines 
should exist, but the intensity of that related to the sec- 
ond array is proportional to both (weak) substrate-array 
and inter-array interaction, so this line is practically not 
observable. The corresponding dispersion curve is plot- 
ted as a thin solid line. On the other hand, the thick solid 
lines are used for dispersion curves corresponding to the 
substrate plasmons and excitations of the first array. 




0.2 0.4 0.6 0.8 
FIG. 7. Frequencies of substrate plasmons and modes of 

QCB with interaction (solid lines) and without interaction 

(dashed lines) as functions of azimuthal angle ip for q — 0.03Q. 



Dispersion curves along the line ABC, (-^ 



2<p _ 
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0, 0.5, 1.0) that includes the triple resonance point B, 
are displayed in Fig. 7 (here q 0.03Q). Here at point 
B, all three modes are strongly coupled and the frequen- 
cies are described by Eq.(38). As in the previous figure, 
dispersion curve corresponding to the second array plas- 
mons is plotted by a thin solid line. 



D. Landau Damping 

Dielectric losses of an isolated substrate are described 
by an imaginary part of its dielectric function e s (22). 
This imaginary part is non-zero within the dashed region 
in Fig. 4 due to appearance of imaginary part of the bare 
polarization operator IIo(k, w) (20) 
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As was mentioned above, the plasmon dielectric losses 
are related to the Landau damping of collective excita- 
tion with momentum k > k* , k*rs ~ 0.1 into electron- 
hole pair. 

The substrate-QCB interaction remarkablly changes 
the conventional picture of substrate plasmon dielectric 
losses. The imaginary part of the dielectric function 
e(k, u>) (18) of the system differs from zero not only 
within the dashed region. New regions of Landau damp- 
ing appear due to QCB - substrate interaction. Indeed, 
QCB serves as a diffraction grid for the substrate like 
in the case considered in Section III. The corresponding 
Umklapp processes initiate Landau damping of surface 
plasmons within the long wave region k < k* where the 
usual Landau damping is absent. The main contribution 
is related to the renormalization term w(m) in Eq.(30) 
due to Umklapp processes along the xi axis (summa- 
tion over rri2 in the expression for function F(m\) is im- 
plied). It is proportional to the fourth power of QCB - 
substrate interaction W 4 . The Umklapp processes along 
both directions xi t 2 contribute also to the renormaliza- 
tion term <p(mi) in Eq.(29). However, they contain addi- 
tional small parameter <p 4 related to inter-array interac- 
tion within QCB and we do not take them into account. 

These new Landau damping regions can be described 
in detail with the help of a "phase diagram" displayed in 
Fig. 8. Coordinate axes in the figure correspond to the 
surface plasmon momenta fci j2 along directions mea- 
sured in units. The QCB period is a = 30 nm. Such 
a choice enables us to realize a rich variety of possible 
damping scenarios. 
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FIG. 8. "Phase diagram" describing possible types of new 
regions of Landau damping. Superlattice constant is chosen 
as a — 30 nm. For an isolated substrate Landau damping 
occurs outside the arc PNTMUR. New regions are bounded 
by curves SIT and WJQU. 

The Landau damping region of the isolated substrate, 
defined by \k\rs > 0.1, is situated outside the circular 
arc PNTMUR. A new damping region inside this arc 
is bounded by two curves SIT and WJQU. The line 
LBJIN (KBQM) corresponds to resonance between 
substrate plasmon and the first (second) array QCB plas- 
mon. The points /, J, T, Q, M, U of intersection be- 
tween the resonant lines, the curves SIT and WJQU, 
and the arc PNTMUR define six rays OI, OJ, OT, OQ, 
OM, and OU and corresponding six angles <p\ = 16°, 
p 2 = 33°, p> 3 = 38°, p> 4 = 69°, ip h = 76°, p 6 = 81°. Each 
pair of adjacent rays (together with coordinate semiaxes) 
bounds a specific structure of new damping region. 

New damping region corresponding to the substrate 
plasmon propagating within the first sector p < ipi is 
separated quite well from the region of initial Landau 
damping. In fact it presents a new damping band. Its 
boundaries are defined by inersection point of the ray 
with a fixed angle <p> with the curves SIT and WJQU. 
The damping amplitude is small because of the small fac- 
tor of order W A mentioned above. When the ray tends 
to the ipi direction, small peak appears near the "blue" 
boundary of the new damping region (precursor of the 
resonance between the substrate plasmon and the first 
array QCB plasmon). 

Within the second sector p>\ < (p < tp 2 , well pro- 
nounced resonant peak appears within the damping 
band. It is described in detail in Fig. 9 for ip = 20°. 
The peak amplitude is of the order of damping ampli- 
tude within the initial damping region (panel a). It has 
a well pronounced Lorentzian form (panel b) placed on 
the wide and low "pedestal" (panel c). 

The next two sectors ip 2 < p < P3 and 933 < ip < ip^, 
do not contain resonant peaks at all. The new damp- 
ing region corresponding to the first one of them is still 
separated from the initial damping region touching it at 
only at the angle ips . For ip 3 < p < p>4, the new damping 
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band is altered by a weak damping tail. The damping 
amplitude in the tail is small due to the same reasons 
mentioned above. It is displayed in the insertion to the 
Fig. 10 for the angle <p = 68°. This angle is close to the 
sector boundary and a precursor of the resonant peak is 
well pronounced. 
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(a) 



this case one deals with a resonance between the sub- 
strate plasmon and the QCB plasmon in the second ar- 
ray. Existence of this resonance is caused by inter-array 
interaction that brings additional small parameter to the 
imaginary part of the dielectric function. As a result, the 
width of the peak is extremely smaller in the second sec- 
tor while surprisingly the peak amplitude has the same 
order of magnitude as in the case of resonance with the 
nearest to substrate (first) array. 

Within the sixth sector p 5 < <p < tp e the damping tail 
does not contains any resonance. It vanishes at <p = ipe 
and within the last seventh sector ip$ < ip < tt/2 the 
damping occurs only within the initial landau damping 
band. 
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FIG. 9. Additional damping region for ip — 20°. a) main 

peak and initial Landau damping region; b) central part of 

the peak zoomed; c) total new damping region and edge of 

the initial damping region 
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0.05 0.06 0.07 0.08 0.09 0.1 0.11 
FIG. 10. Damping tail for ip = 68° and the initial damping 
region. Precursor of the resonant peak is resolved quite well. 

Further increase of the angle p^ < <f < p>5 leads to 
re-appearance of the resonant peak within the tail. In 



The existence of the additional QCB band (tail) of 
Landau damping and appearance of the resonant peak 
within the band (tail) is a bright manifestation of inter- 
play between really 2D surface plasmons and quasi-2D 
QCB plasmons. 



V. CONCLUSION 

In conclusion, we investigated the possibility of spec- 
troscopic studies of the excitation spectrum of quantum 
crossbars, which posesses unique property of dimensional 
crossover both in spatial coordinates and in (q, uS) coordi- 
nates. It follows from our studies that the plasmon exci- 
tations in QCB may be involved in resonance diffraction 
of incident electromagnetic waves and in optical absorp- 
tion in the IR part of spectrum. 

In the case of direct interaction of external electric field 
with QCB, infrared absorption strongly depends on the 
direction of the wave vector q. One can observe dimen- 
sional crossover from ID — > 2D behavior of QCB by 
scanning an incident angle. The crossover manifests it- 
self in the appearance of a set of absorption doublets 
instead of the set of single lines. At special directions, 
one can observe new type of crossover where doublets re- 
place the single lines with changing frequency at a fixed 
q direction. 

Capacitive contact between QCB and semiconductor 
substrate does not destroy the LL character of the long 
wave excitations. However, quite unexpectedly the in- 
teraction between the surface plasmons and plasmon-likc 
excitations of QCB essentially influences the dielectric 
properties of a substrate. First, combined resonances 
manifest themselves in a complicated absorption spec- 
tra. Second, the QCB may be treated as the difraction 
grid for a substrate surface, and Umklapp diffraction pro- 
cesses radically change the plasmon dielectric losses. So 
the surface plasmons are more fragile against interaction 
with superlattice of quantum wires than the LL plasmons 
against interaction with 2D electron gas in a substrate. 
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Dimensional crossover in QCB plays a significant role 
in all the above phenomena. 
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